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Abstract 

We present new and streamlined proofs of various formulae for products and ra- 
tios of characteristic polynomials of random Hermitian matrices that have appeared 
recently in the literature. 
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1 Introduction 



In random matrix theory, unitary ensembles of N x N matrices {H} play a central role 
[in]. Such ensembles are described by a measure da with finite moments |x|'^c?a(x) < oo, 
k = 0,1,2, ■■■ , and the distribution function for the eigenvalues {xj = Xi{H)} of matrices H 
in the ensembles has the form 

rfP«,Ar(x) = -^A{xfda{x) (1.1) 

N 

where da{x) = Y\d,a{xi), A{x) = Y\ ~ ^j) is the Vandermonde determinant for 

i=l N>i>j>l 

the Xj's, and = J ■ ■ ■ J A{xyda{x) is the normalization constant. The special case 
da{x) = e~^^dx is known as the Gaussian Unitary Ensemble (GUE). For symmetric functions 
f{x) = f{xi, ■ ■ ■ , xn) of the Xi's, 

</L = ^ / ■■■ / f{x)A{xrda{x) (1.2) 



N 



denotes the average of / with respect to dPa,N- 

Recently there has been considerable interest in the averages of products and ratios of 



N 



the characteristic polynomials D^lfi, H] = ]^(/i — Xi{H)) of random matrices with respect 

i=l 

to various ensembles. Such averages are used, in particular, in making predictions about the 
moments of the Riemann-zeta function, see fSl ^1 (circular ensembles) and [Sj (unitary 
ensembles). Many other uses are described, for example, in pp, jT2] and [TT] . 
By ()1.2|) . for unitary ensembles, such averages have the form 



Y[-UD^[e„H]/^ ZnJ J n-iiRf 



-A{xfda{x). (1.3) 



In this paper we consider certain explicit determinantal formulae for ()1.3|) - see ()2.6p . ()2.24|) . 
(|2.36p . (j3.3|) . (j3.12j) below. Formula (|2.6p is due to Brezin and Hikami 3j (see also |,17j . 
and when all the fi/s are equal, see [ID]), whereas (|2.24|) . (|2.36|) . (|3.3p and (|3.12|) are due to 
Fyodorov and Strahov fT^ ITT] . The papers ^21 ^] also contain a discussion of the history of 
these formulae. The formulae ()3.3p and ()3.12|) are particularly useful in proving universality 
results for the ratios ()1.3|1 in the Dyson limit as — oo (see |Hj). For a discussion of 
other universality results, particularly the work of Brezin-Hikami and Fyodorov in special 
cases, we again refer the reader to jTTj. The asymptotic analysis in ^I] is based on the 
reformulation of the orthogonal polynomial problem as a Riemann-Hilbert problem by Fokas, 
Its and Kitaev [H]. The Riemann-Hilbert problem is then analyzed asymptotically using the 
non-commutative steepest-descent method introduced by Deift and Zhou |HI, and further 
developed with Venakides in [7] to allow for fully non-linear oscillations, and in [U], [H]. 
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Our goal in this paper is to give new, streamlined proofs of ()2.6|) - ()3.12|) . using only the 
properties of orthogonal polynomials and a minimum of combinatorics. Along the way we 
will also need an integral version of the classical Binet-Cauchy formula due to CAndreief 
dating back to 1883 (see Lemma I^^Tl below) . 

Let TCjiz) = + ■ ■ ■ denote the j'th monic orthogonal polynomial with respect to the 
measure da, 

7Tj{x)7Tk{x)da{x) = CjCkSjk, j, k>0, (1.4) 



where the norming constants c/s are positive. The key observation in our approach is that 
for K = 1 and M = in ((O)) 

{D4fi,H])^ = Mfi) (1.5) 
(see IHj). In our words, the orthogonal polynomial ■nN^jj) with respect to da is also precisely 

N 

the average polynomial ]^(/U — a;j) with respect to dV^^N- Formula p.Sj) appears already in 

i=l 

the work of Heine in the 1880's (see [IH|)- Set 

da^'^^\t)^^^^^^^^^da{t), £,m>0, (1.6) 

(o?a[°'°](t) = da{t)), and let T:f"'\t) denote the jth monic orthogonal polynomial with respect 
to (iat^'™'. With this notation we see immediately from (11.31) . (11.51) that (^tt-^tt-^^') 

is proportional to 7rj^~^'^^'(/ii^) Using a classical determinantal formula of Christoffel (see 
HE]) for 7rK'°l(/i) and a more recent formula of Uvarov [T^ for '7r|^'™'(/i), we are then led 
(see Section 121 Formulae of Christoffel- Uvarov type) to ()2.6|) . ()2.24|) and ()2.36|) in a rather 
straightforward way. Formula ()3.3p appears to have a different character from ()2.6|) . ()2.24|) . 
()2.36|) . and relies on Lemma lTTl mentioned above, which computes the integral of the product 
of two determinants: formula (j3.12p follows (see Section |3| Formulae of two-point function 
type) by combining (|3.3p with (j2.6p and (|2.36p . In [TT] the authors present a variety of 

additional formulae for ( |^m^ )a cases of K and M not covered by ()2.(i|) - ()3.12|) : 

we leave it to the interested reader to verify that the method of this paper can also be used 
to derive these formulae in a straightforward manner. 

Remark 1.1. As is well-known (see e.g., U^), each measure da gives rise to a tridiagonal 
operator 

/ai 6i \ 

hi a2 62 



J = J{da) 



62 03 '■ 

V •••/ 



bi>0 (1.7) 
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with generahzed eigenfunctions given by the orthonormal polynomials 

= c"^vrj(a;), j = 0, 1, ■ ■ ■ , (1.8) 
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i.e. 



hj-ipj-i{x) + ajPj{x) + bjPj+i{x) = xpj{x), j >1 (1.9) 

where bo = 0. Conversely, modulo certain essential self-adjointness issues, da is the spectral 
measure for J in the cyclic subspace generated by J and the vector ei = (1, 0, 0, ■ ■ ■ )"^ (see, 
e.g., |lj). It follows that the transformation of measures 

da da^^^"^^ (1.10) 

leads to the transformation of operators 

J{da) Jida^^^"^^). (1.11) 

For appropriate choices of /ii, ■ ■ ■ , fim and ei, ■ ■ ■ ,ee, such transformations corresponding to 
removing m points from the spectrum of J{da) and inserting i points: in the spectral theory 
literature, such transformations are known as Darboux transformations. The formulae in this 
paper clearly provide formulae for the generalized eigenfunctions p^^'™'(x) of the Darboux- 

transformed operator J{da^^''^^), as well as the matrix entries, a^^'™' and b^^'"^\ in terms of 
the corresponding objects for J{da). Again we leave the details to the reader. Here the 
elementary formulae 



blida) = - + ^Z^ida)Z^^,{da) ^ d 



log/^ (1.12) 

i=0 ^n+lK.U'O-t) 



where dat{x) = e^^da{x), are useful. 

Technical Remark 1.2. Formulae ()2.6|) - ()3.12|) clearly do not make sense for all values of the 
parameters. In all the calculations that follow, we will assume that da has com- 
pact support, support(da)= [—QiQ], say, and that the Hi's and tj's are distinct 
real numbers greater than Q: under these assumptions, da^^''^^{t) becomes, in particular, 
a bona-fide measure, etc. By analytic continuation one sees that the formulae remain true 
for complex values of {fii} and {e^}, as long as they remain distinct. Furthermore, if the /Xj's 
and e/s are distinct, and Im{ej) ^ for all j, then we can let Q ^ (yo and so the formulae 
are true for measures da with unbounded support. Finally we can, for example, let fij fik 
for some j ^ k, which leads to formulae involving derivatives of the vr^'s, etc. 



2 Formulae of Christoffel-Uvarov type 

We use the notations da, ttj, da^^''^\ t^^-'^\ ... of Section^ In addition, in all the calculations 
that follow we assume that da, {^ij}, {cfc} satisfy the conditions described in Technical 
Remark 11.21 above: the natural analytical continuation of the formulae obtained to complex 
values of the parameters, and the limit Q ^ oo,is left to the reader. 

The following result of Christoffel (see |18|) plays a basic role in what follows. 
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Lemma 2.1. Consider the measure da^^^^\t) = Ylif^j -t) dc^it), where i = 1,2,.... Then 

the n^^ monic orthogonal polynomial TTn'^\t) associated with the new measure da^^'^\t) can 
he expressed as follows: 



7r„(/ii) ■ 


■ ■ 7r„,+^(/ii) 






vr„(t) 




7r„(/ii) . 


■ 7r„+£_i(/ii) 


7r„(/i£) . 





(2.1) 



Proof. Set 





7r„(/ii) ■ 






7rn(/i<?) • 






vr„(t) . 


7r„+£(t) 



(2.2) 



We note that qn'^\t) satisfies the condition J P qn''^\t)da{t) = for all j G {0, . . . ,n — 1}. 

is a polynomial of degree at most n. 



Also g|f'°](/i,) = 0, J 
Now observe that 



and so 



(/il - t) . . . (/i£ - t) 



cia[^'0](t) =0, < J < n 



(2.3) 



which means that qn'^^ (t) divided by the product (yUi — t) . . . [fie — t) is proportional to the 

th 



n 



monic orthogonal polynomial 7rn''^\t) associated with the new measure da'^^'^^{t). Now 



gn'^'(t) cannot vanish for any t = ji^j^i > Q, ^ {/ii, ■ ■ ■ Indeed, if gn'°'(/i£+i 

then there exist {a^jf^Q, not all zero, such that p{t) = Yli=o'^i''^n+i(t) vanishes at 
Thus pit) = p(t)/ Y[i=iif^i — ^) is a polynomial of order < n, and as above, pit) is orthogonal 
to P, < j < n, with respect to the measure da'^^~^^''^^{t). Thus p{t) = and hence 
ao = ■ ■ ■ = ae = 0, which is a contradiction. Replacing £ by £ — 1, we conclude that 



0, 



vr„(/ii) ... 7r„+^_i(;Ui) 

TTnifJ'e) ■■■ 7r„+£_i(/i£) 



(2.4) 



Taking the limit t — >■ cx) and noting that the coefficient of the highest degree of vrif '"^ (t) should 
be equal to 1, we find the coefficient of proportionality and establish formula (j2.H) . □ 

Representation ()2.1|) for the monic orthogonal polynomials associated with the measure 
da^^'^^{t) immediately leads to the following result: 
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Corollary 2.2. The product of monic orthogonal polynomials 11^=0 d^fii^^d with 

respect to the different measures is given by the formula 

7r„(/ii) ■ • ■ 7r„+^(/ii) 



i=o 

where A{fi) = Yl ifJ'i - l^j) 

£+l>i>j>l 



A(/i) 



7rn(/i£+i) ■ ■ ■ ■nn+i{l^i+i) 



(2.5) 



We observe that Corollary ()2.2j) gives the identity for the average of products of random 
characteristic polynomials obtained first by Brezin and Hikami 0. 

Theorem 2.3. Let Dn[^,H] he the characteristic polynomial of the Hermitian matrix H. 
The following identity is valid: 



1 



A(/i) 



7riv(/ii) ... 7r7v+L-l(/Wl) 



(2.6) 



\j=l 

where the average is defined by ()1.2|) . 

Proof. To prove formula ()2.6|) we use the representation for the monic orthogonal polynomials 
in the case L = 1 given in (jl.Sp . 



VrAr(Ai) 



N 



(2.7) 



i=l 



Let zfj-'^^ be defined by 



A\x)da^'^'>\x), £=1,2, 



(2.^ 



AT 

where da^^'^^{x) = Y[ da^^''^^{xi). With this notation, we have 

i=l 



llD4^^,,H] 

\j=i 



ylL,0] r^[L,0] y[L~l,0] y[l,0] 



y[L-lfi] y[L-2fi] 



(2.9) 



Equation ()2.7|) implies that rcn ^'^\j^e) can be represented as the ratio ^'^\ where 

7rj^'°'(/i) = vrAr(yu), and = Z^. Thus we obtain 



L-l 



The above equation together with Corollary ()2.2j) proves formula ()2.6p . 



(2.10) 



□ 
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Remark 2.4. Notice (see equations ()2.7|) and ()2.10|) ) that the average of products of charac- 
teristic polynomials can be rewritten as a product of averages. Namely, 



L \ L 



(2.11) 



where (. . .)„[j-i,o] means the average defined by equation ()1.2|) but with respect to the new 
measure da^~^'^^{x), and da{x) = da^^'^^{x). 

The formula of Christoffel (equation ()2.Hl ) enables us to construct the orthogonal poly- 

e 

nomials associated with the measure da^^'^\t) = Y[ if^j ~ t)da{t) in terms of the orthogonal 

i=i 

polynomials associated with the measure da{t). Now we derive a formula due to Uvarov 
[T^ expressing the monic orthogonal polynomials vriJ''™' (t) associated with the measure 

m 

again in terms of the monic orthogonal polynomials 7r„(t) 
associated with the measure da(t). 

Lemma 2.5. Suppose < m < n. The monic orthogonal polynomials 7r|[''™'(t) associated 
with the measure da^'^''^^(t) can be expressed as ratios of determinants, 



hn—m (*^m) 



hJei) 



7r„(t) 



^n-m(ei) • • • ^n~l(ei) 



(2.12) 



hn—mij^ni) ■ ■ ■ ^n— l(^m) 

Here the hk{tj) 's are the Cauchy transformations of the monic orthogonal polynomials iTkit), 



hk{ej] 



Proof. Set 



g[°'-l(t) 



=-/ 


■Kk{t)da{t) 


2m J 


t-ej 


^n— m(^l. 


■ ■ ■ K{^i, 


hn—m i^^m 


) • • • hnij^ra 


'^n—m (t) 


■ ■ ■ 7r„(t) 



(2.13) 



(2.14) 



Now gn'"''(t) is proportional to the n}^ monic orthogonal polynomial Tin''^\t) with respect 
to the measure da^^'"''\t). Indeed, first observe that 

gl?'™Ht) 



[0,m] 



t-ej 



da{t) = 0, j = 1, ■ ■ ■ , m. 



(2.15) 
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Also, for < A; < n, 



nr.i(E«-«) 



E 



ee-t 



(2.16) 



for suitable constants {f3e} and for some polynomial of degree < n — m. But for < A; < n, 



£=1 



t- ei 



da{t) + / 



(2.17) 



The terms in the sum are zero by ()2.15|) and the final integral is zero by the construction 
^^TT^ of g^""' {t) and the fact that deg p{t) < n—m. Thus (ln™'\t) is proportional to Tin'^\t). 
An argument similar to the proof in Lemma f2. II that 



7r„(/ii) ... 7r„+<?_i(/ii) 



7^0, 



(2.18) 



7r„(/i£) ... 7r„+^_i(/i^) 

shows that the denominator in ()2.12|1 does not vanish. Letting t ^ oo in ()2.14|1 . and matching 
leading terms, we prove Lemma f2. 51 □ 

Remark 2.6. In |19j, Uvarov obtains formulae for Tin'^\t) of type ()2.12|) also in the case 
m > n. These formulae can be used to obtain analogues of (|2.24|) and (|2.36|) below in the 
case M > N. 

Remark 2.7. As noted in [T2tlTl]. the Cauchy transformations /ifc(e) of the tt^'s occur explic- 
itly, together with the tt^'s, in the solution of the Fokas-Its-Kitaev Riemann-Hilbert problem 
for orthogonal polynomials jH]. 

Lemma ()2.5j) implies the following analogue of the Christoffel formula for the Cauchy 
transforms of monic orthogonal polynomials. 

Corollary 2.8. Let /i|°'™^(e) be the Cauchy transform of the monic polynomial n^^'^^t) with 
respect to the measure da^^'"^^{t), 



.[0,m] 



Im J t — e 

Let also < m < n. Then hn'"^\e) has a representation similar to that for the monic 
orthogonal polynomials 7r|i''°' (t) ( equation ()2.1|) ), 



h 



n—m y^l I 



h. 



hr 



h„(e) 



(e - e^) • • • ( 



hfi-mi^l) ■ ■ ■ ^n-l(ei) 



hn—mi^^m) ■ ■ ■ l(^m) 



(2.20) 
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Proof. The above representation follows from formula ()2.12|) and from the fact that 



m+l 



{t-e. 



m+l, 



£— n— ■ 

j = l J kytj 



(2.21) 



Indeed we find from formula ()2.12j) that hn'^\e) is the ratio of the determinants. The 
elements of the last row of the determinant in the numerator are the integrals 



T(n-k{t)da{t) 



2m J {t-e){t-em)...{t-ei)' 
Using identity 1)2.211) and noting that the only term 

1 1 



< k < m. 



t-e (e - e^) . . . (e - ei) 
of the sum ()2.2ip contributes to the determinant, ()2.2(Jj) follows. 

Equation ()2.20|) immediately implies the following analogy of ()2.5|) for the h^^'^^ 's. 



(2.22) 
□ 



Corollary 2.9. Let < m < n. Then the product of the Cauchy transforms of monic 
orthogonal polynomials with respect to the measures m can he written as 

a determinant, 



i=0 



Afel 



hri—m (^1 , 



/in(ei) 



hn—m^^m+l) ■ ■ ■ ^n(fr?i+l> 



(2.23) 



Now we derive the identity for the average of the product of inverse random characteristic 
polynomials. 

Theorem 2.10. Suppose 1 < M < N and let 'jn = where Cn is the normality constant 

defined by equation fll.4j) . Then we have the following formula 

hN-Aii^i) ■■■ /ijv-i(ei) 



M 



\j=l 



-1) 



nAf-1 
j=N-M 7j 

2 ± 

Afe) 



(2.24) 



Proof. When M = 1, we use the identity ()2.21|) together with ()2.7|) and the relation (see 
e.g., HE]) 



to obtain 



7„_i = -2'Kin- 
9 



(2.25) 
(2.26) 



We rewrite the average in equation ()2.24|) as follows: 



where 



M 

yIo,m] _ 



,[0,M] y[0,M-l] 



^[O.M-l] y[0,M-2] 



7[0,0] 

'N-M 

y[0,0] 



(2.27) 



(2.28) 



zj^'''^ = Zj\j and da^^'^^{x) = da{x). The following relation can be observed from equations 
and (EIISI): 



ry[0,m] 



(2.29) 



^N-K-l 

Inserting this relation in (j2.27p we find 

M 



M 



Our result (j2.24p immediately follows from the above equation and formula (j2.23p . 
We now repeat the above considerations for the case 



-t) 



(2.30) 
□ 

(2.31) 



The first result is a ChristofFel type formula for the measure ()2.3H) . which is due to Uvarov 

Lemma 2.11. Suppose < m < n. Then the monic orthogonal polynomials 7rlf'™''(t) 's with 
respect to the measure da^''^\t) have the following representation: 



7T. 



[e,m] 



it) 



{t - fie) . . .{t - /ii) 



hn- 


-m(ei) 


■ ■ hn+e{ei) 


K- 


m(^m) 


■ ■ hn+i{em,) 




mil^l) ■ 






-mifJ't) ■ 


■ ■ 'n'n+iifJ'e) 




—m(t) 


■ ■ T^n+eit) 


hn- 


-m(ei) 


■ ■ hn+i{tl) 


K- 


m(fm) 


■ ■ hn+l{^m) 




-m(/il) • 


■ ■ T^n+iifJ'l) 









(2.32) 
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Proof. As in the previous cases we define glf (t) to be the determinant in the numerator of 
dOl- Observe that 

g[f'™](/ii) = ... = g[^'™](/i,)=0 (2.33) 



and that 



ei-t 



(i''^\t)da{t) 



0. 



The next steps are the same as in the proofs of Lemma ()2.1|) and Lemma ()2.5|) . 
Corollary 2.12. 



(2.34) 
□ 



K 



\j=i 



T^N+K-li.lJ'l) 



■KN-m{.IJ'k) ■■■ TTN+K-lifJ'K) 



JO,M] 



A(^) 



^Ar-Af(ei) 



hN(ei) 



(2.35) 



Proof. Identity (j2.35p follows from equations (|2.1(jp and (|2.32|) once we note that equation 
fl2.32|) can be rewritten in a similar manner as equation ()2.5p . □ 

Finally we generalize Theorem ()2.3j) and Theorem ()2.1()j) and obtain a formula for the 
average of ratios of characteristic polynomials. 

Theorem 2.13. Suppose < M < N . Then the average of ratios of characteristic polyno- 
mials of N X N Hermitian matrices H is given by the following formula: 



M(M-1) 

j=N-AI 



A(^)A(6) 



^Af-A/(ei) 

hN-ui^M) 

t^n-m{i^'i) 

T^N-MilJ-K) 



hN+K-l{^l) 
hN+K-li^-u) 

T^N+K-lifJ'x) 



(2.36) 



Proof. Let a^^'^' = a, Zn''^^ = Zn- Then we have 



AK,M] 

'Jl 

7 [0,0] 



ylK,M] y[0,M] 
y[0,M] ym 



(2.37) 
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I.e. 

^'TTT ) '(II'^m) (lln-.%,H]) . (2.38) 

We use Corollary ()2.12|) and Theorem ()2.10j) to obtain formula ()2.36p . □ 

Remark 2.14. Observe that formulae ()2.fi|l . ()2.24|1 do not follow immediately as special cases 
of ()2.36|) : some further algebraic manipulation is required. Similarly, the process of adding 
and removing zeros is clearly reciprocal. More precisely, given ei, ■ ■ ■ , e^, we can construct 
the polynomials TTn'^\t\ c/a'^'^l) associated with the measure da^^'^\t) = {Y[i=i{^i ~ dt by 

^IJ^ : We can then construct vrlf'"' (t; with /x^ = e^, inserting 7r^^^(t; for 

7r„(t) on the right-hand-side of (ED). We should find that 7r|f'°^ (t; rf(a[°'^l)[^'°]) = Unit; da). 
However, again, this relation is not immediately clear, and requires further algebraic manip- 
ulation. 

3 Formulae of two-point function type 

The following integral version of the Binet-Cauchy formula is due to Andreief and plays 
a basic role in our calculations. 

Lemma 3.1. Let {X, dfj,) be a measure space and suppose fi, gj G L'^{X, dfi) for I < i, j < k. 
Then 



X J X 



det{fi{xj))i<ij<kdet{gi{xj))i<ij<kdjj{xi) ■ ■ -d^i^Xk) 



(3.1) 



= A;!det| / fi{x)gj{x)djj,{x) 

\Jx / l<i,j<k 

Proof. Set Cij = fj^fi{x)gj{x)dfi{x). Then 

det{fi{xj))i<ij<kdet{gi{xj))i<ij<kdn{xi) ■ ■ ■ diJ,{xk) 



X J X 



^ sgn(or) sgn(r)c^(i)^(i) ■ ■ ■ c„(^k)r(k) 

CT,TeSk 

^ sgn(cr) ^ sgn(r o cr) 0^(1)^0^7(1) • • • c^(fc)roa(fc) '^^•^^ 

cr T 

gn(r)ci^(i) ■ ■ ■ CkT{k) 



k\ det(c 



ij )l<i,j<k 



as desired. In ()3.2p we used sgn(rocr) = (sgn r) (sgn a) and the fact that c„(i)roa(i) ■ ■ ■ Ca{k)Toa(k) 
cir(i) ■ ■ ■ Cfcr(fc) for all a. □ 
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Theorem 3.2. Let K >1. Then the following identity is valid: 

K 



<K 



where 



and 



I,N+K 



'KN+K{x)T(N+K-l{y) - T^N+K{y)'r^N+K-l{y) 



C 



N,K 



n 



x-y 



N+K-l 2 
l=N 



\2K 



(3.3) 

(3.4) 
(3.5) 



{cn+k-1) 

where q is again the norming constant for tt^ given in p.4p . 

Proof. Let Pj{x) = cJ'^iTj(x), j > 0, denote the orthonormal polynomials with respect to da. 
From ()1.2|1 we obtain 

Adding columns, we see that the Vandermonde determinant A(x, A) has the form 

TToixi) 7ri(xi) ■■■ 7rAr+i^_i(a;i) 



7ro(Ai) 7ri(Ai) 



T^N+K-li^l) 



7ro(Ax) 71i{Xk) ■■■ 7Tn+k-i{^k) 



and similarly for A{x,fi). Here TTj{t) = rc^^'^^t). The determinant A(a;, A) can be evaluated 
by a Lagrange expansion of the form 

7rii(Ai) ■■■ vri^(Ai) 



(3,7) 



E 

0<ii<i2<-<ik<N+K~l 



7rii(AA') ■■■ 7ri^-(Ai^) 



where CTj^ ... = ±1 is an appropriate signature and {(ji, ■ ■ ■ , jat) : < ji < j2 < • • ■ < 
Jn ^ N + K — 1} is the complement of {ii, ■ ■ ■ , ix} in {0, 1, ■ ■ ■ , + — 1}. Multiplying 
(j3.8|) by a similar expansion for A(x,/i), and inserting in (j3.6|) . we obtain a sum of terms of 
the form 



TTj^ixN) 



TX 



da{x) (3.9) 
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which is equal by Lemma ITTl to A^! det (/ TTj/_ (x)??^,, {x)da{x)) ^ 
From this we see that 



<i,k<N 



K 



Y[DN[X„H]DN[fij,H] 



ZnA(X)A(u) ^ 



7rii(Ai) 
7rii(Ax) 



7riK(-^i) 



AT 



'Jfc 



k=l 



mil 



ZnA{x, A)A(x, /i) 



7rii(/ii) ••• vri^(/ii) 



N+K-l 2 
q=N S 



(3.10) 



0<ii<---<ik<N+K-l 



N+K-l 2 
q=N S 



det 



Z7vA(x, A)A(x, yu) \ 



f E 



0<i<N+K-l 



where the last line follows by applying Lemma f3. II to the discrete measure dfi = ^j=o 
But by the Christoffel-Darboux formula 

E/^ N / N '^N+K(y^j)'^N+K~l{fJ'k) — '^N+K{fJ'k)''^N+K-li^j) 
P^{Xj)Pi{Hk) = — 

0<i<iV+A'-l ^ 

which then implies fl3.3|) as Zj^ = N\ Y[f=o^ (see, e.g. [TH]). 

Theorem 3.3. Suppose 1 < K < N . Then the following identity is valid: 

K 



N+K-l 



(3.11) 



□ 



n 



_^^KiK-.y.^K_^ ^ A(e,/.) ^ det(M^,,,^(e.,/.,)),<,,<^ 



where 



/j.Af(e)7rjv-i(/x) - /j-jv_i (e)7r jv(/i) 
e — fi 



(3.12) 



(3.13) 



am 



d again hk{e) = J '^'^^^^^^^^^ is the Cauchy transform of kit) and ^n-i = — 27ri/C^_]^. 
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Observe first that by linearity 



^Ar-M(ei) 
7rAr-Af(/ii) 



■ T^N+L-li^l) 



da{\) 



T!'N-m{^m) 
7rAr-A/(/ii) 

'n'N-uifJ'i) 



7rAr+L-l(Ai) 

7rw+L-i(Ati) 

T^N+L-lifJ'L) 



(3.14) 



Inserting ()2.3(')|1 on the left-hand-side, and using ()2.5|) to re-express the integrand on the 
right-hand-side, we obtain the following result, which is of independent interest. The result 
expresses averages of ratios of characteristic polynomials in terms of averages of products of 
such polynomials. 

Proposition 3.4. Let I < M < N . Then 



da{X) 



X 



M 



(3.15) 



Proof of Theorem\E3 For M = L = K < N , by (prTH|l and dSH, 



A(/.)A(e) 



(-i)KiK-m n;i-_^7, \ nf=i D^ie,, m 



da{X) 



K K 



(3.16) 



[znt) tj) i=i j=i 
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But 



2TTi J Xj - tj Xj - fi^. 



i=l 

K 



"J 3 / \j=2 
K 



2t\% 



i^k 



(3.17) 



as / rfa(Aj)A^7rAr_i(Aj) = / da{Xj)X^j-nN{Xj) = 0for0<£<i^-2<iV-l. Continuing in 
this way, the integral reduces to Hilil/^j ~ ej)M///,Ar(ei, A^fc)- Thus we find 

A(;.)A(6) /nf=i^^[/^.,^]\ A(e,/i) 

det(vy/,Ar+i^(Aj,/ifc)j 



l<i,fc<X 



(3.18) 

and (jTT^ follows. □ 
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